Non-Perturbative Aspects of Fano Resonances in Quantum Dots: An Exact Treatment 
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We consider transport through quantum dots with two tunneling paths. Interference between 
paths gives rise to Fano resonances exhibiting Kondo-like physics. In studying such quantum dots, 
we employ a generalized Anderson model which we argue to be integrable. The exact solution is non- 
perturbative in the tunneling strengths of both paths. By exploiting this integrability, we compute 
the zero temperature linear response conductance of the dot and so obtain reasonable quantitative 
agreement with the experimental measurements reported in Gores et al. PRB 62, 2188 (2000). 
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In recent years advances in nanotechnology have made 
possible the fabrication of single electron transistors 
(SETs). SETs, or colloquially, quantum dots, are char- 
acterized by the remarkable ability to tune, via a gate 
voltage, the number of localized electrons sitting upon 
the dot. By tuning the number of electrons to be nearly 
one, a novel Kondo system is created. This new realiza- 
tion of an old physical paradigm has sparked tremendous 
experimental interest, i.e. ||fi|. 

In the simplest realization of quantum dots there ex- 
ists a single conduction path through the dot. However 
the flexibility inherent in semi-conductor SETs allows for 
more exotic scenarios. Dots may be fabricated such that 
multiple tunneling paths are activated. Generically the 
presence of multiple tunneling paths leads to interference 
effects observable in transport properties. In the pres- 
ence of two tunneling paths, one resonant (i.e. energy 
dependent), one not, Fano resonances can arise [|. Such 
resonances in quantum dots have been observed |QJ^] in 
the form of asymmetric peaks and dips in the linear re- 
sponse conductance as a function of the gate voltage. 

Although this letter will focus upon the observations 
of [Q, Fano resonances are ubiquitous in nanodevices. 
Quantum dots have been embedded in multiple con- 
nected geometries permitting precise delineation of pos- 
sible tunneling paths. Such geometries, in addition, al- 
low threading by Aharonov-Bohm fluxes. The behaviour 
of transport in such devices has been studied in ||. 
Fano resonances have also been observed in STM mea- 
surements of adsorbed magnetic atoms on metallic sub- 
strates. Here interference occurs because of an interplay 
between the Kondo resonance and tunneling into the con- 
tinuum of surface conduction electrons |^,|| . 

The appearance of Fano resonances in quantum dots 
occurs in conjunction with Kondo-like phenomena. Fano 
resonances reported in shows a logarithmic depen- 
dence upon temperature reminiscent of the Kondo effect. 
In addition the authors of 0] observe a sharp depen- 
dence of the Fano resonances upon small magnetic fields. 
Although attributed to a loss of coherent transport in 
the resonant scattering channel, it might also represent 
the destruction of a putative Kondo effect in the dots. 



Observations of Fano resonances in STM tunneling ex- 
periments H are directly related to a Kondo resonance 
arising from the proximity of a magnetic adatom. 

To model Fano resonances we generalize the standard 
two lead Anderson model in the simplest possible way by 
adding a direct lead-lead coupling. In the standard model 
electrons transit from one lead to the other through hop- 
ping on and off the dot. The direct lead-lead tunneling 
provides a competing scattering path, nominally inde- 
pendent of energy. The model Hamiltonian takes the 
form 

it — -^Anderson "b -^lcad— lead tunneling 

/CO 
dxcj IJ (x)d x ci cr (x) 
i=i^,n; er=|,4. '°° 

+V d i(cj cr (0)d cr + h.c.) + e d ^^n a + Un^ni 



+V LR (clj0)c Rr7 (0)+h.c). 



(1) 



Here Ti. encodes the standard Anderson model together 
with an additional term allowing electrons to transit di- 
rectly from one lead to the other. The ci^'s/da's spec- 
ify the lead/dot electrons with n a = d' a d a . U measures 
the Coulomb repulsion on the dot while e^ gives the dot 
single particle energy. Vu are the dot-lead hopping ma- 
trix elements. Experimental realizations of quantum dots 
generically see Vhd ^ Vr<i- Vlr marks the strength of the 
direct transmission channel. The spatial variable x runs 
from -co to co reflecting the 'unfolding' of the leads fl4) . 

In employing this model to describe the transport 
properties of dots such as those studied in Q , we assume 
that only a single level on the dot is relevant to transport. 
This requires the level broadening, Y ~ V 2 , to be consid- 
erably less than U + Ae, where Ae is the level spacing. 
At least for a subset of the data in GJ], this condition is 
met with Y/(U + Ae) ~ .1. We note that in experimen- 
tal measurements on dots with a single tunneling path 
||, Y / (U + Ae) ~ 1/6, yet the Anderson model does an 
excellent job of describing the scaling behaviour of the 
reported finite temperature linear response conductance. 

The model also presumes that the second tunneling 
path is simply due to lead-lead hopping. While it is not 
entirely clear the quantum dots of ¥E are so described, 
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we will show this model more than adequately describes 
the observed phenomena. For dots embedded in a mul- 
tiply connected geometry, i.e. the ambiguity is lifted 
and this term provides a precise description of the second 
tunneling path. Fano resonances in quantum dots were 
described using random matrix theory in H where the 
exact nature of the direct path need not be specified. As- 
pects of the observations in Ref. jl| were described by this 
treatment. However Coulomb interactions were unable to 
be dealt with directly. As we will discuss we believe an 
exact treatment of the non-perturbative physics present 
at finite U is necessary to describe the observations. 

The two lead Anderson model with lead-lead tunnel- 
ing has been studied previously fl(],|ll| . There the model 
was analyzed by expressing all the relevant correlators in 
terms of the dot Greens function (dJd) via a system of 
Dyson equations. The dot correlators are then computed 
via an equations of motion technique j[o| , or a numerical 
renormalization group |p"l| . In our exact treatment of the 
model we find qualitative differences with this approach. 
We believe this is a result of the non-perturbative physics 
inherent in the problem. Although the Dyson equations 
sum up all diagrams, they assume nonetheless the prob- 
lem to be perturbative in the lead-lead coupling, Vlr- 
Our Bethe ansatz solution indicates this to not be the 
case. In particular for U > 0, we do not find a smooth 
Vlr — > limit. Perhaps this is not so unsurprising: we 
similarly have no expectation that the problem is pertur- 
bative in the dot-lead coupling, Vdi- 

Analysis of model: We first examine the particular 
case of VdR = VdL- We will later show how this can be 
generalized. To argue that the model is thus integrable 
we recast the two lead Anderson model into an even/odd 
basis via the transformation, 



c e / = (cl ± cr)/V2. 
With this transformation the Hamiltonian becomes, 
H = H e + H 

/CO 
dxc\ a (x)d x c ea (x) +e d ^n a 
-°° a 
+Un ]ni + r x / 2 (4(0)d CT + h.c.) + V LR c\ a c e<J \ x=0 ; 

U = -«V / dxcl a (x)d x c oa (x) - V L rc 0(t c 00 -\ x =o, (2) 



with T = {V$ L +V$ R ) 1/2 the total dot-lead coupling. Note 
that the dot only couples to the even electrons. 

In changing to an even/odd basis we are still able to 
compute scattering amplitudes of electronic excitations 
off the dot. The even-odd excitations we employ scatter 
off the dot with a pure phase, S e /S . The corresponding 
reflection (R) /transmission (T) probabilities of an elec- 
tronic excitation in the original basis are given by 



Because of the simplicity of the odd sector, S is en- 
ergy independent and given by S = 2tan~ 1 (V / L^). The 
zero temperature linear response conductance is given in 
terms of T by 



G = T = 



(e + q) 2 



where by identifying e = cot(<5 e /2 + taxv' 1 {Vlr)) and 
q = — cot(2 tan^ 1 (Vlr)) we have recast G in a Fano-like 
form. We now turn to the non-trivial computation of 8 e . 
To compute S e we argue H e is solvable via Bethe 



ansatz. To do so we proceed as in 13 for the ordinary 
Anderson model. As a first step we identify an appro- 
priate basis of single particle excitations with momenta 
{kj}. These single particle eigenstates scatter off the dot 
with a bare phase S(k) — — 2tan~ 1 (r(/c — e^) -1 + Vlr)- 
We then proceed to compute the scattering matrices of 
these excitations via computing two particle eigenstates. 
These scattering matrices are identical to that of the ordi- 
nary Anderson model. In particular they satisfy a Yang- 
Baxter relationship. As such multi-particle eigenstates 
can be constructed in a controlled fashion. For a set of N 
particles, their momenta, {kj}, must satisfy the following 
quantization conditions 



N 

n 



^ikj L+i8{kj ) 



A Q — g[kj) + i/2 
X a - g{kj) - i/2 



g{kj) - X a + i/2 



M 

JJ </(*;) "A,,- t/2' 
M \ \ I 



P=l 



Aa — Af 



(4) 



T/R = \(e ld ° =F e ld °)\ 2 /4. 



(3) 



where g(k) = (k - e d - U/2) 2 /2UT. The M auxiliary 
parameters, {A Q }, arise in forming a multiparticle state 
carrying total S z = (N - 2M)/2. The integrability of 
H. ei a new result, leads to a set of quantization conditions 
identical to that of the original Anderson model but for 
one difference: 8(kj) has a different form. This difference 
however is determinative of the physics. 

To compute the dressed scattering phase, S e , of an elec- 
tronic excitation, we employ an argument used by An- 
drei in computing the T = magnetoresistance arising 
from magnetic impurities in a bulk metal |l5f| . The mo- 
mentum, p, of an added electron is determined by the 
quantization conditions of a periodic system of size L via 
p = 2'Kn/L. This momentum has two contributions, one 
coming from the bulk of the system and one from the 
dot, i.e., 

p = 2lin/ L = pbulk + Pimpurity/i- 

The contribution coming from the dot, necessarily scaling 
as 1/L, is to be identified with the scattering phase of the 
excitation off the dot, i.e. S = p- lmp . 

The components of the momenta, Pbulk/Pdotj are 
dressed by the fact that the ground state of the dot-lead 
system is a filled Fermi sea of N interacting electrons 
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(interacting inasmuch as the dot U is finite). In the lan- 
guage of the Bethe ansatz, an N particle ground state 
for e d > -U/2 is formed from N total k's, N - 2M 
of them real, the remaining 2M complex. The 2M 
complex k's are given in terms of M real A a 's with 
each A specifying two k's via k± = x(X) ± iy(X), with 
x(X) = U/2 + e d - VUT(\ + (A 2 + l/4) 1 /2)i/2 ) y(A j = 
\/c7r(-A + (A 2 + l/4) 1 /2)i/2. 

Under the Bethe ansatz the k's are not to be thought of 
as bare momenta of electrons. Rather the Bethe ansatz 
affects a spin-charge separation with the k's associated 
with charge excitations and the A's with spin excitations. 
To compute an electronic scattering phase we must glue 
together contributions coming from the spin and charge 
sectors Q . In adding a spin f electron to the system we 
both add a real k excitation as well as a hole in the set 
of A-excitations. The electronic scattering phase is then 
given by 



5} 



PL 



spin 
im 



p(A). 



The method of computing impurity momenta is discussed 
in detail in Jh|. The impurity momenta are related in 
turn to the impurity densities, Pi mp (k) / a- lm p(X) , of the 
k/X excitations via 



dkP 



charge 



d xPl 



-2lT<T; r 



-imp -Vimpi ^Ai^imp — "" u mip- 

Pimp and (Timp are then governed by the equations, 



p imp (fc) = A(k) + g'(k) / d\ ai (g(k) - A)cr imp (A); 
a imp (A) = A(A) - / dX'a 2 (X' - A)a imp (A') 



— D 



a%(X - g(k))pi mp (k), 



(5) 



where A(A) = — d\5(x(X) + iy{X))/'K and a n (x) — 
2n/ir(n 2 + 4a; 2 ). Q/B mark the 'Fermi surfaces' of the 
seas of k and A excitations while Q is related to the band 
cutoff, D. For the purposes of this paper we are only 
interested in computing the scattering of electrons at the 
Fermi surface. At the Fermi surface, 8\ is given by 



p^ik = B)+ P »^(X = Q). 



Fcrmisurfacc 



The scattering of spin J, excitations can be handled via a 
particle-hole transformation |H. 

We point out that A(A) does not have a smooth 
Vlr — ► limit, a notable difference with the results found 
in Ref. IlQjLlt] . We thus do not expect the problem to per- 



turbative in Vlr- We also emphasize that Pi mp and <7i mp 
encode all degrees of freedom scaling as 1 /L [L is the sys- 
tem size) including corrections to the conduction electron 
density, and not merely those living on the dot. We now 
compute the linear response conductance at T = 0. 



Linear response conductance at H — 0: In pL two 

well developed Fano resonances as a function of the gate 
voltage are reported. The resonances, plotted in Figure 1, 
appear as asymmetric dips. To model these resonances, 
we need to take into account VdR ^ Vdh- To implement 
the even/odd basis change we then use c e / = (V^l / r c l± 
V(IR/lCr) I V^- To assure the even and odd sectors do not 
interact we must allow an additional term to appear in 
the Hamiltonian, SH = (V] L — V] r )(c l cl — c' r cr)/T. 
This term produces lead-lead backscattering. For weak 
asymmetries between Vdh and VdR-, it should not unduly 
affect the physics. 
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FIG. 1. A plot describing the T — linear response conduc- 
tance. The solid curve corresponds to that predicted by the 
Bethe ansatz given the parameters T = .0517, Vtd/Vijd = .75, 
and Q = 0.78, while the circles correspond to experimental 
data reported for a pair of Fano resonances reported in Gores 
et al. [2]. 

We now determine the necessary parameters, U, V, 
VdL/Rt and Vlr entering the model Hamiltonian. From 
the spacing of the peaks together with their widths, the 
ratio Y/U is given by 1/20. To determine the values 
of VLd/VRd and Vlr we use the fact that G for large 
values of the gate voltage tends to its U = value of 
7 46 2 /(l+e 2 ) 2 , with 7 = 4y d 2 L V d 2 ? /r 2 and 9 = V LR /^j. 
Furthermore the value of Vlr determines the depth of the 
dip in G. Using then only two data points we fix these 
later two parameters. 

With these in hand, we manage to reproduce the linear 
response conductance for the entire span of both peaks 
(see Figure 1). We note that the data to which we com- 
pare our theoretical calculations was taken at a low but 
finite temperature. We expect the extremely sharp fea- 
tures seen in our T = computation to be washed out 
at this temperature. 

A feature of these resonances is that G never van- 
ishes. In terms of our computation we believe this to 
be the result of our non-perturbative treatment of fi- 
nite Coulomb interactions. In the free case G(e d ) = 
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7 sin 2 (tan -1 (6 — T/cd) + tan _1 (0)) and G vanishes for 
some ed- Similarly the expression for G arising out of the 
Dyson equations [lO|JTl[ | always vanishes for some value 
of €d, one reason we suspect that the Dyson equations 
do not adequately capture the physics at finite U, Vlr- 
We also point out that the resonances occur for values 
of the gate voltage placing the dot in its mixed valence 
regime (n<j < 1) and not the Kondo regime (rid ~ !)• 
Generically, our solution predicts that the linear response 
conductance in the Kondo regime will be relatively struc- 
tureless. 
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FIG. 2. A set of Fano resonances for differing values of F. 
These curves are computed using Vlr = 0.4. 

Dependence of width of Fano resonances upon 

T: In Q, Fano resonances were studied as a function 
of the total dot-lead coupling strength, T, where it was 
observed that the width of Fano resonances exhibit a non- 
monotonic dependence upon T. (In a dot with a single 
tunneling path, the width of a resonance merely increases 
with r.) Together with this non-monotonicity, the overall 
shape and amount of asymmetry in the Fano resonances 
was observed to be sensitive to the strength of T. 

We can reproduce this array of behaviour. Plotted in 
Figure 2 is the linear response conductance for a set of 
differing T's. For T small, a Fano resonance appears as 
a sharply peaked bipolar structure. As T is increased, as 
naively expected, the bipolar peak broadens. However at 
some critical value of T ~ (.3 — A)U, the bipolar reso- 
nance is replaced by a narrow unipolar one. With further 
increases in T, this resonance proceeds to broaden out. 

Linear response conductance at H ^ 0: The be- 
haviour of Fano resonances in magnetic fields was also 
studied in 0. It was found that the resonances exhib- 
ited a marked response to extremely small magnetic fields 
(giiH/T ~ 1CP 2 ). In particular they demonstrated that 
upon application of H, a small bipolar Fano resonance 
is transformed into a much larger unipolar structure (see 
inset to Figure 3). We are able to reproduce this phe- 



nomena (see main body of Figure 3). For H = 0, a small 
bipolar resonance in G is present. With the introduction 
of a small field, a large unipolar peak is superimposed 
over the bipolar structure. As this calculation is done at 
T = 0, finite T should lead the two structures to merge 
leaving a reasonable representation of the experimental 
data. 

The strength of H necessary to produce the unipolar 
peak is on the order of a putative Kondo temperature, 
Tk, which at the symmetric point of a single channel dot 
is estimated by, T k - v / l7rexp(-C//8r). This might 
suggest that in applying H, a resonance (or lack thereof) 
due to the Kondo effect is destroyed. 
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FIG. 3. The response of a bipolar Fano resonance to the 
application of small magnetic fields, H. To compute these 
curves we employ Y/U = 1/20 and Vlr = —3.6. Inset: Ob- 
served response of a Fano resonance to small applied fields. 

To summarize we have employed a generalized Ander- 
son Hamiltonian to describe observations of Fano res- 
onances in quantum dots. We have argued that this 
Hamiltonian is integrable and sketched how this integra- 
bility can be exploited to compute the T = linear re- 
sponse conductance. Using this model, we are able to de- 
scribe a number of observed features of Fano resonances 
presented in Our exact solution of the model sug- 
gests that the physics underlying the resonances is non- 
perturbative in the presence of finite Coulomb repulsion 
on the dot. 
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